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Abstract 

It has been shown recently that deformations of some integrable field theories in (1 -|- 1)- 
dimensions possess an infinite number of charges that are asymptotically conserved in 
the scattering of soliton like solutions. Such charges are not conserved in time and 
they do vary considerably during the scattering process, however they all return in the 
remote future (after the scattering) to the values they had in the remote past (before the 
scattering). Such non-linear phenomenon was named quasi-integrability, and it seems to 
be related to special properties of the solutions under a space-time parity transformation. 
In this paper we investigate, analytically and numerically, such phenomenon in the 
context of deformations of the integrable Bullough-Dodd model. We hnd that a special 
class of two-soliton like solutions of such deformed theories do present an infinite number 
of asymptotically conserved charges. 



1 Introduction 


The objective of the present paper is to investigate, analytically and nnmerically, the con¬ 
cept of qnasi-integrability, hrst proposed in [1], in the context of deformations of the exactly 
integrable Bnllough-Dodd model [2]. The motivation of onr study is to try to shed light on 
the mechanisms responsible for such interesting non-linear phenomenon which has a large 
potential for applications in many areas of physics, mathematics and non-linear sciences in 
general. 

As it is well known, solitons in (l-l-l)-dimensions are solutions of non-linear held equations 
which travel with constant speed without dispersion and dissipation, and when they scatter 
through each other they keep their forms, energies, etc, with the only effect being a shift in 
their positions relative to the ones they would have if the scattering have not occurred. Such 
extraordinary behavior is credited to the fact that the solitons appear in the so-called exactly 
integrable field theories in (1 -|- l)-dimensions, that possess an infinite number of exactly 
conserved charges. Therefore, the only way for the scattering process to preserve the values of 
such an inhnity of charges, is for the solitons to come out of it exactly as they have entered it. 
In addition, in most of such theories the strength of the interaction of the solitons is inversely 
related to the coupling constant, i.e. the solitons are weakly coupled in the strong regime 
and vice-versa. Such behavior and the large amount of symmetries (conserved charges) make 
the solitons the natural candidates for the normal modes of the theory in the strong coupling 
regime, opening the way for the development of many non-perturbative techniques in the 
study of non-linear phenomena. The drawback of such approach is that exactly integrable 
soliton theories are rare, and few of them really describe phenomena in the real world. 

The observation put forward in [1] is that many theories which are not integrable present 
solutions that behave very similarly to solitons, i.e. such soliton like solution scatter through 
each other without distorting them very much. It was shown in [T] in the context of deforma¬ 
tions of the sine-Gordon model, and then in other theories mm, that such quasi-integrable 
theories possess and infinite number of charges that are asymptotically conserved. By that 
one means that during the scattering of two soliton like solutions such charges do vary in time 
(and quite a lot sometimes) but they all return in the remote future (after the scattering) to 
the values they had in the remote past (before the scattering). Since in a scattering process 
what matters are the asymptotic states, such theories are effectively integrable, and that is 
why they were named quasi-integrable. 

The mechanisms behind such non-linear phenomenon are not well understood yet. All 
the examples studied so far are deformations of exactly integrable field theories. The zero- 
curvature condition or Lax-Zakharov-Shabat equation [B] of the integrable theory becomes 
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anomalous when applied to the deformed theory, and so the Lax potentials fail to become 
flat connections when the equations of motion hold true. Despite those facts, techniques of 
integrable held theories can be adapted and applied to construct an inhnite number of charges 
Q^^') which present an anomalous conservation law 




( 1 . 1 ) 


The anomalies have some interesting properties. They vanish exactly when evaluated in 
one-soliton type solutions, and also vanish for two-soliton type solutions when the two solitons 
are well separated. The anomalies are only non-vanishing when the soliton like solutions are 
close together and interact with each other. For some special classes of soliton like solutions 
the anomalies have a further striking property. They have a mirror type symmetry in 
the sense that the charges have the same values when reflected around a particular value of 
time ^A, which depends upon the parameters of the solution. In other words, one hnds that 
(t^ = (“^)’ "where i = t — t^. So, the charges are not only asymptotically conserved, 

i.e. (oo ) = (— 00 ), but are symmetric with respect to a given value of the time. The 

only explanation found so far for such behavior of the charges, is that those special soliton 
like solutions transform in a special way under a space-time parity transformation, where the 
point in space-time around which space and time are reversed depend upon the parameters 
of the solution under consideration. The proof of the connection between parity and mirror 
symmetry of the charges involves an interplay of the Lorentz transformations and internal 
transformations in the Kac-Moody algebra underlying the anomalous Lax equation. We do 
not believe however that such parity property is one of the causes of the quasi-integrability, 
but it seems to be present whenever such phenomenon occurs. 


In this paper we investigate the concept of quasi-integrability in the context of deformations 
of the Bullough-Dodd model |2] involving a complex scalar held in (1 -|- l)-dimensions with 
Lagrangian given by 

^ = ( 1 . 2 ) 

and the potentials being 

(1.3) 

with e being a real deformation parameter, such that the Bullough-Dodd model is recovered 
in the case £ = 0. Because of some particularities of the vacuum solutions of such models, 
as explained in section the physically interesting deformed theories exist only when the 
parameter e: is restricted to rational values. 


We construct the anomalous zero-curvature condition for the theories (1.2) with the Lax 
potentials taking values on the twisted sl{2) Kac-Moody algebra . The charges 
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satisfying (1.1) are obtained by the so-called abelianization procedure [3 El [HI Uni E] where 
the Lax potentials are gauge transformed into an inhnite abelian sub-algebra oi A 2 . In fact, 
due to the anomaly of the zero-curvature only one component of the Lax potentials can be 
rotated into such sub-algebra, leading therefore to the anomalous conservation (1.1). The Lax 
potentials do not transform as vectors under the (1 -|- l)-dimensional Lorentz transformations. 
However, the grading operator of the Kac-Moody algebra generates a one-dimensional 
subgroup isomorphic to the Lorentz group, and we show that the Lax potentials are vectors 
under the combined action of those two groups. That fact allows us to show that the anomalies 
in (1.1) do vanish when evaluated on the one-soliton solutions of the theories (1.2). In addition, 
we show that some special two-soliton solutions lead to the existence of a space-time parity 
transformation P such that the complex scalar held (p, when evaluated on them, transforms 
as P {(p) = p*. For such two-soliton solutions we show that the real part of the charges 
satisfy a mirror symmetry, as described above, and are therefore asymptotically conserved. 
The imaginary part of the charges however, are not asymptotically conserved. 

We also implement a perturbative method to construct solutions of the deformed theory 
(1.2) as a power series in the deformation parameter e, as p = p^ + e pi + p 2 + ■ ■ -i such 
that pq is an exact solution of the integrable Bullough-Dodd model. We then split the helds 
into their real and imaginary parts and then into their even and odd parts under the parity 
transformation P, i.e. Pn^^’^ = ^ (1 -P) ^ with p^ = p^ + By starting with an 

exact solution pq of the Bullough-Dodd model that satishes P ( 920 ) = or P {‘■Po') — P ’0 
and P (^Pq^ = —Po, we show that the pair of helds , Pi~^ satisfy a pair of linear non- 

homogeneous equations, and the pair of helds [pi~ , Pi^^ satisfy a pair of linear homogeneous 
equations. Therefore, it is always possible to choose solutions where [pi~ , Pi^^ = 0, and 
so the hrst order held satishes, P {pi) = p*i- Once that is chosen, one can show that the 
same structure repeats at the second order in the expansion and one can choose solutions 
such that P {p 2 ) = P 2 - By repeating such procedure order by order we show that the theories 
(1.2) always contain solutions with the property P {p) = p*, and so the charges evaluated 
on them present the mirror symmetry described above. So, the dynamics of the deformed 
theories (1.2) favors the “good” solutions, in the sense that it is not possible to have solutions 
satisfying the pure opposite behavior under the parity, i.e. P {p) = —p*. That is an interesting 
interplay between the dynamics and the parity that deserves further study. For instance, the 
production of “bad” modes could be energetically disfavored and emission of radiation could 
be suppressed. 

We also perform numerical simulations, based on the fourth order Runge-Kutta method, to 
study the scattering of two soliton like solutions of (1.2). We performed simulations for various 
rational values of the deformation parameter £, and found that in all cases the predictions 
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of the analytical calculations were confirmed, i.e. the real part of charges do satisfy the 
mirror symmetry when we use as a seed for the simulations solutions of the Bullough-Dodd 
model that have the right parity property, i.e. P {(f) = (f*. So, the evolution of the fields 
under the deformed equations of motion seem not to destroy the parity property of the initial 
configuration, again indicating that the dynamics seem to favor the “good” modes as observed 
in the analytical perturbative expansion mentioned above. The mirror symmetry of the charges 


were checked in the simulations by evaluating the first non-trivial anomaly (see ( 1 . 1 )) as 
well as its integrated version dt' (3^^^ = (t) — ( — 00 ). All simulations show 

that the real part of 7 ^^^ is symmetric under reflection around a given value of time close to 
t = 0, and so leading to the mirror symmetry for the real part of the charge and then 
for its asymptotic conservation. The imaginary part of 7 ^®^ is not symmetric under reflection 
and does not lead to the asymptotic conservation of the imaginary part of the charge 

The paper is organized as follows: section [^discusses the properties of the vacuum solutions 


of the theories ( 1 . 2 ) and their implications on the possible physically interesting deformations 


of the Bullough-Dodd model. In section we present the construction of the quasi-conserved 
charges using techniques of integrable field theories based on the anomalous zero-curvature or 
Lax-Zakharov-Shabat equation. The interplay between the Lorentz and parity transformations 
leading to the mirror symmetry of the charges is dicussed in section]^ and section|^implements 


the perturbative method to construct solutions of ( 1 . 2 ) as power series in the deformation 


parameter e and discusses the connection between dynamics and parity. The Hirota’s one- 
soliton and two-soliton exact solutions of the integrable Bullough-Dodd model are given in 
section The numerical simulations are presented in section and our conclusions are given 
in section]^ The appendix [A| gives some basic results about the twisted si (2) Kac-Moody 
algebra A 2 used in the text. 


2 The deformed Bullough-Dodd models 

We shall consider models of a complex scalar field (p in (1 -|- l)-dimensions with Lagrangian 


given by (1.2) and the potentials being given by (1.3). The Euler-Lagrange equation following 


from ( 1 . 2 ) is 


dfcp-dlcp + e^= 0 ( 2 . 1 ) 

where we have taken the speed of light to be unity. The real and imaginary parts of the 


equation ( 2 . 1 ) are given by 

- dl^pR + cos ((p/) - cos ((2 -)- e) (p/) = 0 

d^(pi-dl(pi+ e'^^ sin {(fj)+ sin {{2 +e) (fi) = 0 


( 2 . 2 ) 
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where we have denoted + i The Hamiltonian associated to (1.2) is conserved and 

complex, and denoting the Hamiltonian density as Ti = Hr + iHi, we get 

Hr = ^ {dt<PRf + {d^p>Rf - + Vr 


Hi — dtp>R dtp>i + dxP>R dx<pi + Vj 
with V = Vr + iVj, and 


(2.3) 


U,^ = 

cos (pi) + 

^ p-(2+£)¥3H 

2 + £ 

U/ = 

sin {pi) — 

^ p-(2+£)v7i? 

2 + £ 


(2.4) 


Note that the densities (2.3) are not positive dehnite, and so we can not really define vacuum 
configurations as those having minimum energies. However, in order for the space integrals 
of the densities (2.3) to be conserved in time, one needs the momenta to vanish at spatial 
infinity, and so one needs the fields to be constants there. Therefore, from the equations of 

.at SI 

dVi 


motion (2.2) one gets that such constant configurations are extrema of the potentials, i.e. 
OVr 


dipR dipi 

dVj dVR 


dp}R dipi 

which implies 

glS+elv’i? gQg ^ gQg ('('2 -I- g'j 


= cos {(pi) — e cos ((2 + e) ipi) = 0 

= sin (ipi) + sin ((2 + e) ipi) = 0 


(2.5) 


g(3+£) Vr gjj^ ^ — sin ((2 + e) p>i) (2.6) 


Squaring both equations and adding them up one concludes that (pR = 0. Using that fact and 
manipulating ( |2.6 ), by multiplying and adding them up, one concludes that sin ((3 + e) pi) = 
0 and cos ((3 + e) pi) = 1. Therefore, the extrema of the potentials are 

27rn' 


{pR , Pi) = (^0 , 


3 + £ 


with n integer 


(2.7) 


Note that such extrema are not maxima or minima of the potentials. They all correspond to 
saddle points. 

We now come to a very interesting property of such models. For static configurations we 
have that the quantities 

1 r. 


Sr — —^ {^xPr)'^ — (dxPi)'^ + Ur ; 


£i — —dxPR dxPi + Vi 


( 2 . 8 ) 


are constant in a: as a consequence of the static equations of motion (2.2), i.e. ^ = 0. 

Indeed, those quantities correspond to the static Hamiltonian densities (2.3) with the space 
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coordinate x replaced by an imaginary time ix. So, the quantities (2.8) are conserved in 


the time ix for a mechanical problem of a particle moving on a two dimensional space with 


coordinates (pR and ipi. Therefore, for a given static solution of the model (1.2) that goes. 


at spatial inhnity, to vacua conhgurations (2.7), one concludes that the potentials Vr and Vi 


(and so the complex potential V) are bound to have the same values at both ends of spatial 


inhnity. If one denotes by and n_ the integers labeling the vacua (2.7) at x —)■ oo and 
X —)■ —cx), respectively, then from (1.3) one gets that 


i 2 TT —'n,_ ) 

e (3+=) = 1 


(2.9) 


For the Bullough-Dodd model where e = 0, one gets that n+ — n_ has to be a multiple of 
3. Indeed, the static one-soliton solutions of that model do not “tunnel” between consecutive 
vacua as x varies from —oo to +oo, like in the sine-Gordon model, but jumps 2 vacua and 
ends in the third one. 


For the deformed Bullough-Dodd models where £ 7 ^ 0, the only way of satisfying the 


condition (2.9) for any real value of £ is to have n+ = n_, i.e. for a static solution the vacua 


are the same at both ends of spatial inhnity. If one wants non-trivial one soliton solution 
with non vanishing topological charge, then the deformation parameter e has to be taken to 
be a rational number. That is a quite striking and interesting restriction on the ways the 
Bullough-Dodd model can be deformed. If one takes 


P 

2 + e=- 
Q 


with p and q integers 


( 2 . 10 ) 


then the equations of motion ( 2 . 1 ) takes the form 


aU - sill, + e”!" - e-’"^ ^ 0 


( 2 . 11 ) 


where we have redehned the helds as (p = q(j), and the space-time coordinates as a; = ^x, 
and t = ^Jqt. 


3 The quasi-conserved charges 

We now use techniques of exactly integrable held theories in (1 -|- 1) dimensions to construct 
quasi-conserved charges for the non-integrable theories we are considering. We introduce the 
connection (Lax potentials) 

71+= (-21/ + m) 61 + ^Fi ; A_ = - d_ cp Fq (3.1) 

dtp 

where we have used light-cone coordinates as x± = and so d± = dtF dx- In addition, 
the operators bn and Fn are generators of the twisted loop (Kac-Moody) algebra A 2 dehned 
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in appendix [A} Even though we will be interested in potentials of the form (1.3), in the 


connection (3.1) we shall assume that E is a generic potential that depends upon the scalar 


held (p but not on its complex conjugate p*. In addition, we shall assume that the potential 
does not involve complex parameters in its dehnition, such that its complex conjugate can be 
obtained just by the replacement p ^ p*, i.e. 


V* = V{p^ p*) 


(3.2) 


The reasons for assuming that property will become clear when we discuss below the anomalies 
in the conservation of the charges. 


One can check that the curvature of the connection (3.1) is given by 

d+A_ -d-A+ + [A+, A_] = - (^d+d-p + Fq - d-pXFi 

with 


d^V dV 

^ ~ o —2 o- 2V + m 

op^ op 


The coefficient of Fq in (3.3) corresponds to the equation of motion of the theory 


dV 

a+a_(p + —= 0 

op 


(3.3) 


(3.4) 


(3.5) 


and, when it holds true the vanishing of the curvature depends upon the vanishing of the 
anomalous term X. Note that by shifting the potential asV —)■ E + one observes that X 
vanishes only when E is a linear combination of the exponential terms and e'^. So, the 


curvature (3.3) vanishes for the Bullough-Dodd potential, corresponding to (1.3) for e = 0, 


or then for the potential of the Liouville model, E ~ e‘^. The fact that the Liouville model 

/2'\ 

admits a zero curvature representation in terms of the twisted Kac-Moody algebra A 2 is 
perhaps not know in the literature. 

In order to construct the charges we use the so-called abelianization procedure u\mm 
[Tn] or Drinfeld-Sokolov reduction m- We perform a gauge transformation of the deformed 


connection (3.1) as 


A^^ a^ = g A^g ^ - d^g g (3.6) 

with g being an exponentiation of the positive grade generators introduced in the appendix 

0 
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^n — Cn Fn 


(3.7) 


Splitting things according to the grading operator (A.5), we have that the a_ component of 


the transformed potential (3.6), becomes 

00 

E (n) 

; 


n =—1 


j (A 
a, a_ 


(n) 

= na_ 


(3.8) 
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The components of it are 
= 6-1 

a® = -[6-1, J-i]-a_(pFo (3.9) 

= -[6-1, J-2] + 2 [[6-1, 


aL^ — — [6-1, J^n+i ] + • • • 

The crucial algebraic property used in the following calculation is the fact that 6 -i is a semi- 
simple element of the algebra ^ = Ag p in the sense that Q is split into the kernel and image 
of its adjoint action and they have an empty intersection, i.e. 

^ = Ker -I- Im [ 6 -i,Ker] = 0 Im=[ 6 -i,^] (3.10) 

From the appendix one observes that the elements of Ker are bQn±i, and the elements of 
Im are F^. Since Fn+i is in the image of the adjoint action of 6 -i, none of its components 
commute with 6 -i . Therefore one can recursively chooses the parameters C^+i inside the Fn+i 
to kill the component of in the image, i.e. in the direction of After that procedure is 
done the connection O- takes the form 

U- = 6-1 -|- 6 i -|- 65 -|- ... (3-11) 


The hrst two non-trivial components are 


(3.12) 


and 


A-(5) _ ^ 


^ {d-^) - Y i^-^) (9-^) 
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5 {d-ipf + ^ (al(p) {d_^f + ^ (a!(p) {d_^) 


d^ip{d-ip) 


(3.13) 


Once the parameters (n are chosen to rotate O- into the abelian subalgebra (kernel) generated 
by bQn±i, there is nothing we can do about the a+ component of the transformed connection 
(3.6). We only know it will have positive grade components only, since g is generated by the 
positive grade F^s. Therefore, it is of the form 




N=1 


n=l 


N = Ql±l-, 


/ = 0,1,2,... (3.14) 



In the case of integrable field theories where the quantity X, dehned in (3.4), vanishes, it is 
possible to show that by using the equations of motion the image component of a+ vanishes. 
In our case, the use of the equations of motion (3.5) can show that all the quantities are 
linear in X, given in (3.4). 


When the equations of motion (3.5) hold true the transformed curvature reads (see (3.3)) 
cl+a_ — cl_a+ + [a+ , a_ ] = —d_(p X g Fi g~^ (3.15) 

We now write 
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F,g-^= + 


N=1 


71 = 1 


and the hrst two components of the hrst term on the r.h.s. of (3.16) are 
= 0 


(3.16) 


(3.17) 


= 


4 {8-<pf - 4 (dlifi) - 2 (d-ifi) (atifi) - 2 (atf) 

Since a_ does not have components in the direction of the F„’s, it follows that the commutator 
in (3.15) does not produce terms in the directions of the ^at’s. Therefore, one has that 


which in terms of the space-time coordinates x and t becomes 




Dehning the charges as 

we get 

dt 




dx 


= piN) 


(3^^^ = - r dxd.ifXa^^^ 

Z J—oo 


(3.18) 

(3.19) 

(3.20) 

(3.21) 


' —OO 

XX) 


where we have assumed boundary conditions such that a”’''"'’'' {x —)■ cxd) — {x —>■ —oo) = 0. 
We call (3^^) the anomaly of the charge Q^^'). An useful quantity in our numerical simulations 
is what we call the integrated anomaly 7*^^^ dehned as 


^X) = qX) {t) _ qX) (_ oo ) = l f dt' r dxd.^Xa^^^ 

2 , j—00 J—00 


(3.22) 


From (3.17) one notes that the anomaly f3^^') vanishes and so the charge Q^^) is conserved. 
It corresponds in fact to a linear combination of the energy and momentum. The hrst non 
trivial anomaly corresponds to = 5 and the expression for and X are given in (3.17) 
and (3.4) respectively. 
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4 The role of Lorentz and parity transformations 


Consider the (1 + l)-dimensional Lorentz transformation 

X — vt 


A 


x± —>■ x± 


or 


X —>■ 


\/r^ 


t —y 


t — vx 
\/l — v‘^ 


(4.1) 


where a is the rapidity and v the velocity, i.e. v = tanha. The Lax potentials (3.1) do not 


transform as vectors nnder snch Lorentz boost. However, consider the automorphism of the 
loop algebra A 2 

S(T) = (4.2) 

where d is the grading operator dehned in (A.5). It then follows that the Lax operators ( |3.1 ) 
transform, under the composed transformation, as vectors, i.e. 


n (A±) = 


± 


with 


H = AS 


(4.3) 


Therefore, the curvature (3.3) is invariant under such combined transformation, and so is the 


anomalous term d-cpX Fi. In order to see how the anomalies of the charges (3.20) transform 


under H, we have to inspect the properties of the quantities introduced in (3.16) 


Note that the term Fq, appearing on the r.h.s. of the second equation of (3.9), trans¬ 
forms under H as 

n {d-(p Fo) = d-(p Fo (4.4) 


As explained below (3.9), we choose the parameter ^1 in Fi such that the term [6_i, Fi] 
cancels the term d^ipF^, and so such two terms have to transform under the same rule under 
H. Since H (&-i) = e“ 6_i, it follows that 


^ {Fi) = Fi 


(4.5) 


Indeed, one hnds that the cancelation implies that one should choose Ci = —d-(p, and so 


H (Fi) = e"Fi, and H ((^i) = e "Ci. Now, using (4.5) one observes that the last three 


terms on the r.h.s. of the third equation in (3.9) gets multiplied by e " under the action 


of H. Therefore, in order to cancel the Fi component on the r.h.s. of that equation one 
needs H ([&_i, F 2 ]) = [&-i, F 2 ], and so one has to have H {F 2 ) = F 2 - Continuing such 

reasoning order by order, one concludes that all Fn’s have to be invariant under H, and so the 


group element performing the gauge transformation (3.6) is also invariant, i.e. 


0.[g) = g 


(4.6) 


Therefore, similarly to A„, the transformed connection a„ also behaves as vector under H, i.e. 


n {a±) = a± 


(4.7) 
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According to the way the -7>i’s are chosen in (3.9) to cancel the components of a_ in the 


direction of the F^s, it follows that the parameters (nS are functions of the x_-derivatives of 


the scalar held ip. But from (4.6) one concludes that 


f!(C..)=A(C„) = e-"“C„ 


(4.8) 


Therefore, each term in (n must contain n x_-derivatives of the held (p. From (4.6) one 
observes that fl {g Fi g~^) = e°‘ g Fi g~^, and so every term on the r.h.s. of (3.16) have to get 
mulplied by e“ under the action of ff. Therefore, since {Bn) = b^, it follows that 


a (a‘">) = A (a<")) 




(4.9) 


From the dehnition of in (3.16), it is clear that it is a function of the C„’s, and so a 


function of the ^--derivatives of the held ip. Therefore, from (4.9) one concludes that each 


term in must contain (N — 1) ^--derivatives of p>. Indeed, from (3.17) one observes that 
contains four ^--derivatives of p>. 

We are now in a position to draw some conclusions about the anomalies of the charges 


dehned in (3.20). Consider a solution of the equations of motion (3.5) which is a traveling 
wave, i.e. (p = (p{x — vt). One can then make a Lorentz transformation and go to the rest 
frame of such solution where it becomes static, i.e. x-dependent only. Clearly the charges 
evaluated on such solution must be time independent, and so its anomaly dehned 


in (3.21), must vanish. But from (3.21) and (4.9) it follows that 


n dt) = A dt 


= e 


-Na 


13^^^ dt 


(4.10) 


Therefore, dt, and so dQ^^\ is a tensor under the (1 -|- l)-dimensional Lorentz group. 
Consequently, if dQ^^'^ vanishes on the rest frame of the solution, it should vanish in all 
Lorentz frames. One then concludes that the charges are exactly conserved for traveling 


wave solutions (like one-soliton solutions) of (3.5). In fact, such conclusion applies for any 


functional of the scalar held ip and its derivatives, which is a tensor under the Lorentz group. 

The one-solitons we treat in this paper are localized solutions in the sense that the held 
ip have non-vanishing space-time derivatives only in a small region of space. Therefore, the 


integrand in the dehnition (3.21) of the anomaly (3^^\ is non-vanishing only in such a small 


region of space, i.e., it gets exponentially suppressed outside such region. In addition, the one- 
solitons interact with each other by short range interactions. So, a two-soliton solution for 
the case when the two one-solitons are far apart should be just the superposition of the one- 
soliton solutions. Therefore, the anomaly evaluated on a two-soliton solution when the 
one-solitons are far apart should vanish, because it reduces to the sum of the anomalies of the 
two one-solitons, which by the argument above vanish. Consequently one should expect the 
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anomaly to be non-vanishing only when the solitons are close together and interacting. 
That is in fact what we observe in onr nnmerical simnlations described in section m. We 
do not have yet a good nnderstanding of non-linear dynamics governing the behavior of the 
charges and anomalies when the solitons interact with each other. That is a crncial issne to be 
nnderstood and is at the heart of onr working dehnition of the concept of qnasi-integrability. 
What is clear so far is that special properties of the solntions nnder a space-time parity 
transformation play an important role in all that. It is not clear however if snch properties 
are the canses or conseqnences of the qnasi-integrability. Let ns explain how it works. 


Let ns then consider a (two-soliton like) solntion of the eqnations of motion (3.5), and a 
space-time parity transformation 


P : 


{x,t) ^ (-X, -t) 


X = X 


Xa 


t = t — t A 


(4.11) 


where the valnes of xa and tA depend npon the parameters of that particnlar solntion. There 
are two important classes of two-soliton solntions according to the way they behave nnder 
the parity transformation. The hrst one is that where the two-soliton solntion is invariant 


nnder the parity, i.e. P {(f) = </?. From the argnments below (4.9) we conclnded that each 


term in mnst contain (N — 1) x_-derivatives of cp. From (3.14) we have that the integer 


N is of the form N = 6 / ± 1, with I any integer. Therefore, each term in contains 
an even number of x_-derivatives of (p, and so it is invariant under the parity. Therefore, 
P {d-^p X = —d-ip X Consequently 


r di r dxd.^Xa^^^ = 0 
J-io J-XQ 


(4.12) 


where Xq and to cire any chosen values of x and t respectively, i.e. we are integrating on 


a rectangle with center in (xa , ^a) (see (4.11)). Therefore, the charges evaluates on such 


two-soliton solutions satisfy the mirror like symmetry 


o'"’ (l„) = (-io) (4.13) 

for any to- The asymptotic conservation of the charges is therefore a particular case of such 
mirror symmetry, corresponding to the case where to —?• oo. 

The physically important two-soliton solutions however, are not invariant under the parity 
transformation. The complex scalar held ip, evaluated on such two-soliton solutions, satisfy 
the property 

P(V)=V>‘ (4.14) 


We are assuming that the potentials V of our theories satisfy the property (3.2), i.e., 
the potential depends only on ip and not on its complex conjugated (p*, and the complex 
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conjugated of V is obtained just by the replacement ip ^ (p*. Therefore the anomaly X, 
dehned in (3.4), also satisfy the same property, i.e.. 


X* =X{ip^ ip*) 


(4.15) 


In addition, the Lax potentials dehned in (3.1), do not involve complex parameters, and 
according to the appendix the structure constants of the loop algebra in the basis 
and Fn, are all real. Consequently, the complex conjugate of the charges and of 

the quantities a'^^\ dehned in (3.20) and (3.16) respectively, are also obtained just by the 
replacement ip ^ ip*, i.e. 


(gW)* = QW(v.^v^*) 






{ip ip*) 


(4.16) 


As we have argued above, each term in the quantity contains an even number of X-- 
derivatives of ip, and so for the two-soliton solutions satisfying (4.14), one gets that 

P{d_ipXa^^^) = -d_ip*X* (4.17) 


Consequently 


j . dt J dx (d-ip X + d-ip* X* ) = 0 


(4.18) 


where, like in (4.12), we are integrating on a rectangle with center in (xa , ^a) (see (4.11)). 
Therefore, for the two-soliton solutions satisfying (4.14), the real part of the charges 
satisfy the mirror like symmetry 


(o'"’ + (y"')*) ko) = (<3'~> + {q'AI (-‘») 

for any to. Consequently, in the limit to —t oo, we get that the real parts of the charges are 
asymptotically conserved, i.e., they have the same values before and after the scattering of 
the solitons. 


5 The interplay between parity and the perturbative 
expansion 

Let us consider potentials V that are perturbations of the Bullough-Dodd potential, in the 
sense that they depend upon a parameter e such that they become the Bullough-Dodd po¬ 


tential for £ = 0. The potential (1.3) is an example of it. We shall expand the solutions in a 
power series in the parameter e as 


ip = ipo + eipi + e ip 2 + ... 


(5.1) 
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Therefore the potential depends explicitly upon e and also implicitly through (p, and so we 
have the expansion 


dV 

dip 


dV 

dip 

2 


£ 


|£=0 

d^V 


d^V d^V dip 


dedip dip"^ de 


(5.2) 


£•=0 


d^V dip d^V d^ ip d^V fdip' 

de‘^dip^ dedip‘^ de dip‘^ de'^ dip^ \de 


+ 


0(.3) 


£=0 


From the equation of motion (3.5) we then get the equations for the components of the held 
ip in the expansion (5.1), as {d‘^ = df — 


with 


d^ip. 


h = 


aVo + e‘^°= 0 

d^V 


dip'^ 

d^V 

dedip 


e=0 fn 


i = 1,2,3... 


|£:=0 


/2 2 


d^V 
de"^ dip 


| e =0 +2 


d^V 

dedip'^ 


, ^ , 2 
e=0 + W—T e=0 

dip^ 


(5.3) 

(5.4) 

(5.5) 


and so on. 


Let us split the helds into the even and odd parts under the parity transformation (4.11), 
and into their real and imaginary parts as well 




± — 


1 


-pi±P); 




R — 


1 


= 2(1 + ^*) ; 




I — 




(5.6) 


with C being the complex conjugation operation. By splitting the zero order equation ( |5.3 ) 
into its even and odd parts under the parity one gets 


^ + 2 




= 0 


(5.7) 


Note therefore that one should not expect non-trivial solutions for the case ip^ = 0, since ip^ 
would have to assume some very special (vacuum) constant values. On the other hand the 
case (po = 0 can lead to non-trivial solutions. Following (4.14) we shall therefore consider 
solutions of the pure Bullogh-Dodd equation (5.3) which satishes 

p (<do) = ; i-e. 


p h?) 


^0 




(5.8) 


Note that /i given in (5.5) is a function of the order zero held ipQ only. But since we are 
assuming that the potential satishes the property (3.2), it follows from (5.8) that the action 
of the parity P on fi is the same as the action of the complex conjugation C. Therefore 

t (1 ± P) 1 (1 ± C) A = ft ; 1 (1 ± P) i (1 T C) A = 0 


(5.9) 
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But since |e=o is also a function of (^o only, it follows from the same reasoning that 


1 1 d'^V 


'5W , ' 

^ 1^=0 


dip^ 


1 1 d'^V 


Therefore splitting the equation (5.4) for /i into its even and odd parts under P and its real 


and imaginary components one gets 

'd‘^V , \ 


+ 

dV{- + 

dVi’^ + 

avf’“ + 


5(^2 

'd^V 

9(^2 

'd^V 

9(^2 

'd^V 


|£ = 0 


| e :=0 


| e :=0 


|£:=0 


-R,+ I • 
(Pi +t 


I- 




R- , • 
(Pl’ +l 


' d^V 

9(^2 

' d^V 

dip^ 

kW 




|£=0 


| e :=0 


| e :=0 


|£:=0 




R,- 


R- 




= ft 


= -ifi 


= 0 


= 0 


(5.10) 


So, the pair of helds satisfy a pair of linear non-homogeneous equations, and 

the pair of helds satisfy a pair of linear homogeneous equations. In addition, 

the two pairs of equations are decoupled. Therefore, [(pt~ i = 0 is a solutions of such 

equations, but the pair ((pf’"'', p>i~^ can never vanish. If one has a given solution (^(pf , ipt) of 
the equations above, then the conhguration {(pi '~^, , ip\^ — [(pi~ , is also 

a solution. Therefore, given a solution one can always make its real part even under P, and 
its imaginary part odd under P, i.e. one can always choose the hrst order solution to satisfy 


P (V7i) = (pi 


(5.11) 


The quantity /2 given in (5.5) is a function of Pq and pi only. If the potential V satisfy the 


property (3.2) then it follows, by the same arguments used above, that 
1 


(l±/’) t ( 1 ±C) f, = f, , 


± . 


t ( 1 ±P) i (ItC) f, = 0 


(5.12) 


Consequently the pairs of helds , pt ) and (p 2 ’ , P 2 ~^^, satisfy equations identical to 


(5.10) with /i replaced by /2. Therefore, by same arguments as above, one can always choose 


the second order solution to satisfy 

P{P2) = PI (5.13) 

Continuing with such process order by order, one concludes that is always possible to choose 


a solution, as long as the perturbative series converge, that satishes the property (4.14) which 
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was used to prove the mirror symmetry property (4.19) satisfied by the real part of the charges 

The conclusion is that the deformations of the Bulough-Dodd model we are considering 
seems to possess a subclass of solutions that present an infinite number of real charges which 
are quasi-conserved. By quasi-conserved we mean charges satisfying the mirror symmetry 
(4.19). In the case of two-soliton like solutions such properties imply that the infinity of real 
charges preserve the same values, after the scattering of the solitons, that they had prior 
the scattering, even though during the scattering process itself they may vary considerably. 
It is that sub-sector of the model, consisting of solutions satisfying (4.14), that we call a 
quasi-integrable theory. 


6 The Hirota’s solutions of the Bullough-Dood model 


We now show that the pure Bullough-Dood model possesses soliton solutions satisfying the 


property (|5.8|). We shall use the Hirota’s method to construct such solutions, and so we 

To 


introduce the r-functions as 


ipo = In 


Tl 


( 6 . 1 ) 


One can check that if such r-functions satisfy the Hirota’s equations 

Tod+d-To - d+To d-To = rf - Tq 
Tid+d_Ti - d+Ti d_Ti = To Tl - 


( 6 . 2 ) 


then the zero order field cpo, given in (6.1), satisfy the Bullough-Dood equation (5.3). 


6.1 One-soliton solutions 


The one-soliton solutions of the pure Bullough-Dood model correspond to the following solu¬ 
tions of the Hirota’s equations (6.2) 


To = 1 — 4 a e -I- a e 


2 „2r 


Tl = (l-f- a 


(6.3) 


with 

T = a/s (^Xj^ — —^ (6.4) 

where z and a are complex parameters. The one-soliton solutions leading to the quasi- 
integrable sector of the deformed theories we are interested in, correspond to the cases where 
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z is real. Then we parameterize it as z = e and define v = tanha, where v is the soliton 

-^/3rE(°) 

velocity and so a is the rapidity. We now write a = e®^e and define ae^ = e'^, with 


Therefore, we get from (6.3) 


W = ^ [x — vt — 


tq cosh IT — 2 


(6.5) 


( 6 . 6 ) 


Ti cosh W + 1 

Note that if the phase ^ vanishes, then ^ vanishes whenever cosh IT = 2, and that corresponds 
to a singularity in the solution for ipQ. Therefore, we shall restrict to the cases where ^ 7 ^ 0. 
We have that 

cosh IT = cos^ cosh ^ (x — vt — + i sin^ sinh 7 (x — vt — 


Therefore, under the space-time parity transformation 


(x — x^^^^ — (x — 


one gets that 


cosh IT —>■ cosh IT* 


and so 


t —y —t 


-t (fo 


(6.7) 

( 6 . 8 ) 

(6.9) 


Therefore, the one-soliton solutions ( 6 . 6 ) satisfy the property (5.8). By the arguments of 


section such one-soliton solution of the pure Bullough-Dodd model can serve as a seed 
to construct, by a perturbative approach, one-soliton solutions of the deformed theory that 


satisfy the property (4.14). 


However, by the arguments presented below (|4.10 ), for any traveling wave solution, and 


so for any one-soliton solution, the charges given in (3.20), are not only quasi-conserved 
but exactly conserved. 


6.2 Two soliton solution 


The solutions on the two-soliton sector of the pure Bullough-Dood model are obtained by 


solving (6.2) by the Hirota’s method, and the result is 


To = 1 — 4 oi e ^ — 4 02 e ^-1-af e ^-f 02 e 
2 zf — zf Z 2 + 2 Z 2 


2 ^ 2 r 2 


+ 801 O 2 


ri+ra 


{Zi + Z2Y {zf + Z1Z2 + Z2) 


, „2 „ (^1 - ^2? {zj -Z 1 Z 2 + zl) ^2ri+r2 

4t Qj-i do o ^ r-. ^ 

{zi + Z2) (zf + Z1Z2 + x|) 
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4ai a 


2 - ^ 2 ? {zf -Z1Z2 + zj) gri+ 2 r 2 

{zi + ^2)^ {zi + Zi Z2 + Z2) 


-\- CL-i d 


2 2 (^1 ^2) (^1 Zi Z2 + Z2) 2ri+2r2 


1 “2 


( 2:1 + 2 : 2 )"^ (zf + Zi Z 2 + Z 2 ) 

Ti = 1 + 2 ai + 2 02 + 02 

zf+ 4 zfzl + 4 


+ 4 oi 02 


Ji+ra 


(^1 + ^2) (^1 + Zi Z2 + Z2) 


o „2 „ (^1 - ^ 2 )^ ^2 + zi) ^ 2 ri+r 2 

^ til CiO O ^ ^ ^ 

(zi + ^2) {zi + Z1Z2 + Z2) 

{zi - ^2)^ {zj -Z1Z2 + zj) gri+2r2 

(zi + ^ 2 )^ (^1 + ^2 + Z 2 ) 


2 oi O 2 


2 2 (^1 Z 2 ) {zf ^ 12:2 + 2 ^ 2 ) ^ 2 ri+ 2 r 2 


Ct-^ ^2 


(Zi + ^2)^ (^1 + Zi Z2 + Z2) 


where 


hi = a/3 


ZiX+ - 


X- 


i = 1,2 


( 6 . 10 ) 


( 6 . 11 ) 


and where Zi and Oj, i = 1 , 2 , are complex parameters. As in the case of the one-soliton 
solutions, we shall deal with the two-soliton solutions where Zi are real, and parameterize the 
solutions with six real parameters Vi, and xf'\ i = 1,2, as 


( 0 ) 


= e-“* ; 


-73 

= V^~'’i 


t;j = tanhQ;j (6-12) 

Similarly to what we have done in the case of one-soliton let us dehne Oj e'"® = with 

i^x — Vi t — 


Wi = Vs 


+ i& 


'1-vi 

In addition we introduce the real quantities A, cq and ci as 


i = 1,2 


(6.13) 


^A(v,) ^ Vl - ^ 2 )^ (zf -Z1Z2 + Zj) 
(Zl + Z2)^ (zf + Z1Z2 + zi) 

2 z\ — zi zi + 2 zi 


Co (Vi) = 


Cl [Vi) = 


(zi - Z2f {zi - 2:1 ^2 + zi) 
zf+ 4 zizi + zj 
{zi - ^2)^ {zi - ^2 + zi) 


(6.14) 
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Therefore, the Hirota’s r-functions become 

n = l + 2e^^ +2e^^ + + e^i+'^2+^ 


8co-4e^^-4e^^+e^^+^^+^ 
4 Cl + 2 + 2 + e^i+'^2+^' 


We then introduce the new space-time coordinates as 


X+ 




Vit — x\^^) (^X — V2 t — 




1-V^ 


+ A 


X 


- ^ ^ 


(^X — Vi t — Xi^^"j (^X — V2 t — X2^^'j 


vi 


(6.15) 


\/^-^i ^ 

Note that as long as vi ^ V 2 , x± ^^^e independent coordinates. In fact, one can check if vi = V 2 


the two-soliton solution (6.10) reduces to a one-soliton solution (6.3) (taking z\ = Z 2 = z and 


oi -|- 02 = a) . In addition we denote 

w+ = Wl + W2 + A = x+ + ^ (6 + 6) ; IT-= 1+1-1+2 = x- + ^ (6-6) ( 6 .I 6 ) 

Therefore, by factoring out from tq and ri, the term 2 one gets that 

tq cosh 1++ -I- 4 Co — 8 cosh cosh + e~^ cosh W- 


Ti cosh W+ -I- 2 Cl -I- 4 e 2 ^ cosh cosh + e“^ cosh l+_ 
Note that under the space-time parity transformation 

P : 


(6.17) 


one gets that 
and consequently 


(x+, X-) (-X+, -X-) 

P{W±) = -l+± 




Tl 


t; 


(6.18) 

(6.19) 

( 6 . 20 ) 


Therefore, from (6.1) one observes that such two-soliton solutions of the pure Bullough-Dood 


model satisfy the property (5.8). By the arguments of section]^ such two-soliton solution can 
serve as a seed to construct, by a perturbative approach, two-soliton solutions of the deformed 


theory that satisfy the property (4.14). 
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7 The numerical simulations 


We performed various simulations of the proposed deformed Bullough-Dodd model to check 
the ideas behind the concept of quasi-integrability, specially the behavior of the anomalies 
(3.21) and (3.22) during the scattering of two solitons. To do so we had hrst to generate static 
solutions to the deformed equations of motion. Then we stitched two of them together to 
perform a kink-kink colision and analysed the results. 


In all our numerical work the time evolution was simulated by the fourth order Runge- 
Kutta method. After several lattice sizes were tested we found that a lattice extending from 
—25 to -|-25 was enough to obtain reliable results. We also found that a grid spacing of 0.01 
in space and of 0.003 in time discretization provided enough accuracy. Our initial condition 
for the kink-kink scattering is not exact, thus radiation was generated. This radiation was 
absorbed at the boundaries through a viscous force acting from x = —25 to a: = —22 and 
from X = +22 to a: = -|-25. Such a boundary condition effectively simulates an infinite grid. 


7.1 Generating the initial condition 


To obtain static solutions of (2.1) we used a minimization method that consists in solving 

g-(2 + 6)sP ^ Q_ 


dr+ - -f e 


(7.1) 


This equation is an evolution in the fake time r. The static solutions of (7.1) and those of 


(2.1) are the same, but only the former is dissipative. The dissipative nature and the absence 


of sources means that any initial condition will lead to a static solution in the infinite-time 
limit. We used the exact solutions of the Bullough-Dodd model as seeds to this method, and 
labeled the results with the original parameters - i.e. seeding the kink obtained with ^ = 0.9 


from the original model (e = 0) into (7.1) results in the deformed kink with ^ = 0.9. 

Once we stitch two kinks together we can use two facts to generate the complete initial 
condition: first, the two kinks (let’s call them the left kink and the right kink) are far from 
each other and can be treated as free kinks; second, each independent kinks have Lorentz 
invariace so 

(Pi = (p,(x-ut) ^ (7.2) 


dx V dt 

where ipi denotes any of the two kinks in the absence of the other. We use ipi as the left kink 


in (7.2) from the grid’s left edge to position x = 0, and as the right kink for the rest of the 


grid. This also introduces some perturbation with respect to the exact two kink solution, but 
one that is small for small colision velocities. 
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7.2 Kink-kink interaction 


Once we have generated our initial conditions we simulated the colision between two kinks, 
and calculated the first non-trivial anomaly and integrated anomaly for N = 5, 


given in (3.21) and (3.22) respectively. Using (3.4), (3.17), (3.21) and (3.22), and for the 


potential (1.3) one gets that (with the choice m = 0 in (3.4)) 
s (£■ + 3) 


/?(5) = I r dx 4:{d-ipf dlif - 4:(dlipY - 2d-ipdlip - 2dt^ 

2 (2 -|- £) J—oo L ^ ' 


(7.3) 


and 


^(5) ^ + 3) 


2 (2 + £) 7 -, 


r / \ 2 

/ dxA{d-ipfd^_ip-A(dlip) - 2d-ipdlip - 2dt^ 

-OO J —oo L ^ ^ 


(7.4) 


Note that e = —3 is special because the potential (1.3) vanishes in that case 


We have performed a large number of simulations of the scattering of two kinks for various 
values of the deformation parameter £ and for the phases i = 1,2, which appear in the 


exact soliton solutions (see (6.5) and (6.16)). Note that such phases are parameters of the 


exact soliton solutions but we do not know how the deformed solitons depend upon them for 
£ 7 ^ 0. However, as explained above, we have constructed the static kink numerically (with 
£ 7 ^ 0) from an exact kink (with e = 0) for a given choice of the phase ^i. So, we still use the 
phases to label the deformed kink. 

The results for the scattering of two kinks for various choices of the parameters £ and 
i = 1,2, were qualitatively the same, and we show in the hgures Biii and 1^ some 
representative cases. On those hgures we show on the top plots, the real and imaginary parts 


of the anomaly given in (7.3), and on the bottom plots the real and imaginary parts of 


the integrated anomaly given in (7.4) 


Note that, in the cases shown in hgures [l]j^ the real part of the integrated anomaly 


7 


( 5 ) 


is symmetric under rehection around a point close to t = 0. From (3.22) one has 


7^^^ (o) - 7^^^ (-O) = (o) - (-O) = ^ / / dxd-ifXa^^'^ (7.5) 

According to the numerical simulations, the real part of ~ 7^^^ ’^^^i^Ties, and 


so one gets a numerical conhrmation of (4.19), for W = 5 at least. Remember from (4.11) 
that i = t — and for the simulations shown in hgures [l](^ one has close to zero. 

The case £ = — 1, shown in[^ has an asymmetry however, in the oscillations of the anomaly, 
which is perhaps due to numerical errors. Therefore, our numerical results conhrm the ex¬ 


istence of a mirror symmetry given in (4.19), and obtained from the analytical calculations 
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based on the parity properties of the solutions. In addition, such results imply that real part of 
the charge is asymptotically conserved, i.e. ReQ^^^ {t = —oo) = Re {t = oo). The plots of 
the imaginary part of the integrated anomaly 7 *-^^ show that the imaginary part of the charge 
is not asymptotically conserved, and so it is in agreement with our analytical calculations 
that lead to the asymptotic conservation of the real part of the charge only. 


8 Conclusions 


We have considered deformations of the exactly integrable Bullough-Dodd model dehned by 


the Lagrangians (1.2) and potentials (1.3), involving a deformation parameter e. We have 


shown that such deformed theories possess sectors where the soliton-like solutions present 
properties very similar to solitons in exactly integrable held theories. They possess an inh- 
nite set of quantities which are exactly conserved in time for one-soliton type solutions, i.e. 
solutions traveling with a constant speed, and are asymptotically conserved for two-soliton 
type solutions. By asymptotically conserved we mean quantities which do vary in time during 
the scattering process of two one-solitons but that return, in the distant future (after the 
collision), to the values they had in the distant past (before the collision). Since what matters 
in a scattering process are the asymptotic states, such behavior is effectively what we observe 
in scattering of solitons in integrable held theories, and that is why we call such theories 
quasi-integrable. The mechanism behind the asymptotic conservation of such an inhnity of 
charges is not well understood yet, but in all examples where it has been observed so far 
[MIlllS], the two-soliton type solutions present special properties under a space-time parity 
transformation. In the case of the deformations considered in this paper, the complex scalar 
held ip of the theory, when evaluated on the two-soliton solutions, is mapped into its complex 


conjugate under the parity transformation (see (4.14)). It is worth mentioning that the point 


around which space and time are reversed, under the parity transformation, depend upon the 
parameters of the solution under consideration. 

The inhnity of asymptotically conserved charges was constructed using techniques of in¬ 
tegrable held theories based on an anomalous zero-curvature (Lax) equation, where the Lax 

/o\ 

potentials live on the twisted sZ(2) Kac-Moody algebra A 2 . We have used the properties of 
the Lax potentials under Lorentz transformation and internal A 2 ^ transformations to prove 


that two-soliton type solution satisfying the property (4.14) under the parity transformation. 


present an inhnity of asymptotically conserved charges. In addition, we have shown that by 


starting with a two-soliton solution of the integrable Bullough-Dodd model, satisfying (4.14), 


one can always construct by a power series expansion on the deformation parameter e, a so¬ 


lution of the deformed model that also satishes the parity property (4.14) and leading to the 
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asymptotic conservation of the charges. So, the dynamics of the deformed model seems to 
favor the qnasi-integrable sector of the theory. 


The nnmerical simnlations have confirmed the predictions of onr analytical calcnlations. 
The nnmerical code was divided in two main parts: the hrst part uses an exact static one- 
soliton solution of the integrable Bullough-Dodd model as a seed to constructed a static 
one-soliton of the deformed model under a relaxation procedure. The second part of the code 
takes two of such deformed one-solitons far apart and stitches them at the middle point, and 
then performs the scattering of them, absorbing the radiation at the edges of the grid. The 
anomalies of the charge’s non-conservation are evaluated to verify their asymptotic conser¬ 
vation. It is interesting to note that the initial conhgurations of the simulations are built 


from exact solutions of the Bullough-Dodd model that satisfy the parity property (4.14). The 


simulations could well drive them into solutions of the deformed model that do not satisfy 


(4.14). However, in all our simulations the real part of the charges were observed to be asymp¬ 


totically conserved indicating that they do satisfy the property (4.14). That corroborates the 


conclusion of our analytical perturbative calculation, mentioned above, that the dynamics of 


the deformed model favors the solutions satisfying (4.14). 


We do not believe that the parity property (4.14) is causing the phenomena we call quasi- 


integrability, but it seems to be present whenever such non-linear phenomena manifest it¬ 
self. We need further investigations to fully understand the mechanisms behind the quasi- 
integrability. That is a quite interesting phenomena that may have several applications in 
many areas of non-linear sciences. 
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( 2 ) 

A The twisted loop algebra A 2 ^ 


We consider here the twisted loop algebra A 2 ; which corresponds to the twisted Kac-Moody 

(2i') 

algebra A 2 with vanishing central term. For onr applications it is convenient to work with a 
realization of such algebra based on SU{3) matrices, and the following basis 




Fen+j = X"Fj J = 0,1, 2,3,4, 5; neZ (A.l) 


where A is the so-called spectral parameter, and 
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UJ = gF2i+l)7r/6 


j = 0,1,2,3,4,5 


(A.2) 


(A.3) 


The commutation relations are given by (independently of the choice of the u given in 


(A.3)) 


[61, 6_i] = 0 


(A.4) 


[bi,F,] 

— Ffc+i 

fc = 0 , 1 , 2 , 3,4 ; 

[b-i,F^] 

= Fi.i 

/ = 1 , 2 , 3 , 4 , 5 ; 

[FcFi] = 

1 

to 

0 - 

1 

[Fi, F2] = -F3 

[Fo,F2] = 

-F 2 

[Fi,F3] = -F4 


[hi,F^] = XF, 
[h.i,Fo]=X-^F, 


24 


[F 2 , F 4 ] — —XFq 
[F 2 , F 5 ] = —2 A 61 




























[ Aq , F3 ] — 0 

[Ai,F 4 ] 

= 2 Xb_i 

[ A3 , F4 ] — 2 A 61 — A Fi 

[Ao,F 4 ] = F4 

[Fi,F,] 

= AFo 

[ A3, F5 ] = —A F2 

[Fo,F5] = -2A6_i + F5 

[F 2 ,F,] 

= —2 A b_i — F5 

[F,,F,] = -XF, 


The sub indices of the generators (A.l) correspond to the grades given by the grading operator 
d = Fq + QX^; [d, Gm] = mgrn 

with 


(A.5) 


6^6n = {A"Fo} 

^6n+3 = 


Gen+i = {X^h,X^F,} 
^6n+4 = {A"'F4} 


gen+2 = {A"F2} (A.6) 

6^6n+5 = {A"+'6-l, A^Fg} 
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Figure 1: Scattering of two kinks for the choice of parameters e = 1, .^i = 0.9 and .^2 = 0.9. 
On the top plots it is shown the real and imaginary parts of the anomaly given in (7.3), 
and on the bottom plots are shown the real and imaginary parts of the integrated anomaly 
^(5)^ given in (7.4). 












Figure 2: Scattering of two kinks for the choice of parameters e = 1, .^i = 0.9 and .^2 = 1-4. 
On the top plots it is shown the real and imaginary parts of the anomaly given in (7.3), 
and on the bottom plots are shown the real and imaginary parts of the integrated anomaly 
^(5)^ given in (7.4). 
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Figure 3: Scattering of two kinks for the choice of parameters £ = 1, .^i = 0.9 and .^2 = 2.0. 
On the top plots it is shown the real and imaginary parts of the anomaly given in (7.3), 
and on the bottom plots are shown the real and imaginary parts of the integrated anomaly 
^(5)^ given in (7.4). 
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Figure 4: Scattering of two kinks for the choice of parameters e = 2, = 0.9 and .^2 = 2.0. 

On the top plots it is shown the real and imaginary parts of the anomaly given in (7.3), 
and on the bottom plots are shown the real and imaginary parts of the integrated anomaly 
^(5)^ given in (7.4). 
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Figure 5: Scattering of two kinks for the choice of parameters e = —1, = 0.9 and ^2 = 2.0. 

On the top plots it is shown the real and imaginary parts of the anomaly given in (7.3), 
and on the bottom plots are shown the real and imaginary parts of the integrated anomaly 
^(5)^ given in (7.4). 
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